JOURNAL OF RESEARCH of the National Bureau of Standards— B. Mathematics and Mathematical Physics 

Vol. 64B, No. 3, July-September 1960 

Electric Polarizability of a Short Right Circular 
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(April 29, 1960) 

A method similar to that employed by Smythe [l] 2 for calculating the capacitance of 
a freely charged short right circular conducting cylinder is used to calculate the electric polar- 
izability tensor in the principal axis system for such a cylinder. Calculations to an accuracy 
of approximately five significant figures are carried out for cylinders with radius to half-length 
ratios of M, }4, 1, 2, and 4. The results are applicable to the design of artificial dielectrics. 

1. Introduction 

A solid conducting object in the form of a short, that is noninfinite in length, right circular 
cylinder makes a suitable element for use in the construction of an artificial dielectric. For 
such an object, therefore, it becomes important to know the electric polarizability tensor 
a tj which relates the induced dipole moment p t to the inducing field Ej according to the equation 

Vi=^E h (1) 

with the cylinder situated in free space. Once a tj is known, the effective dielectric constant of 
a spatial array of identical and identically oriented cylinders can be calculated according to 
the methods of Kaprielian [2]. 

A centered coordinate system in which the 2-axis coincides with the axis of rotational sym- 
metry obviously constitutes a system of principal axes for a short right circular cylinder. 
The components a xx and a yy are clearly equal and are denoted a tt where t stands for "trans- 
verse"; similarly, a zz is denoted a xx where I stands for "longitudinal." The cylinder cross section, 
with nomenclature, is illustrated in figure 1. This article describes a method which has been 
used successfully for calculating a a and a tt for cylinders with radius to half-length (a/b) ratios 
of %j %, 1, 2, and 4. Rationalized M.K.S. units are used throughout. 

The problem which presents itself is that of finding, in the external space, a Laplacian 
function which has no tangential gradient at the surface of the cylinder and which reduces to 
the potential of a uniform field at infinity. With such a potential, there is associated physically 
a surface charge density on the cylinder proportional to the local normal gradient and a field- 
free condition in the interior. It is the presence of this density which gives rise to the induced 
dipole moment. Since there is no known coordinate system in which a cylinder of finite length 
with closed plane ends forms a coordinate surface and in which Laplace's equation is separable, 
the present problem must be regarded as intractable from the point of view of conventional 
methods and another method involving an arbitrarily good approximation to the surface charge 
density must be invoked. 

As a mathematical abstraction, let the cylinder be regarded merely as a geometrical con- 
struct to which the surface charge distribution just mentioned is firmly affixed. If that portion 
of the potential which corresponds to the uniform (applied) field is now subtracted, the remain- 
ing potential will be that due to the charge alone and will be found to correspond to a localized 
external field with dipole-like appearance and a perfectly uniform internal field which is equal 



* This article is based upon a portion of a thesis which has been accepted by the faculty of the California Institute of Technology in partial 
fulfillment of the requirements for the degree of Doctor of Philosophy. (Invited paper.) 
i University of California, Riverside. 
2 Figures in brackets indicate the literature references at the end of this paper. 
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Figure 1. Cross section of the cylinder with nomenclature. 

and opposite to the applied field. The problem may therefore be restated in terms of finding, 
on a finite cylindrical surface with closed plane ends, a charge distribution which, when acting 
alone, generates a uniform field in the interior. The technique employed in the present article 
makes use of this viewpoint and consists in setting up, on the surface, a charge density function 
which is completely determined by a finite number of parameters, then solving for the values of 
these parameters so as to obtain maximal uniformity of the charge field or, alternatively, 
maximal cancellation of the applied field, within. This method is similar to that employed by 
Smythe for calculating the capacitance of a freely charged cylinder [1]. The research reported 
here was performed under the direction of Prof. Smythe, to whom the author is indebted 
for many valuable discussions. The cooperation of the personnel of the California Institute 
of Technology Computing Center is also gratefully acknowledged. 

A second article, now in preparation, applies this method to the problem of determining 
the magnetic polarizability tensor under the assumption of negligible field penetration. If 
both tensors are known and if the wavelength is long compared with element spacing, it becomes 
possible to calculate the wave propagation constant of artifical media whose elements are cyl- 
inders of the type considered here. 

2. Functions for Expressing the Charge Densities 

A sj^stem of generalized orthogonal polynomials and weighting functions has been con- 
structed for the purpose of expressing the charge densities on the side and ends of the cylinder. 
These polynomials involve an argument v, which represents either z/b or p/a as required, and a 
weighting function (1— u 2 )\ The parameter v is given the value of minus one-third in order 
that the charge densities may be asymptotic to a constant times Z _1/3 as Z— K), where I is the 
distance from the edge of the cylinder. The product of the polynomial of order m and the 
weighting function is denoted by the symbol "ty m and is defined as follows: 



lW7,$>,^)= 



(_ l)m2->r(m+Y-r-f) 

= r( 7 +f)r(m+i-H 
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(2) 



If the parameter 7 is given the value of one-half, this function becomes suitable for use on the 
side ol the cylinder; if unity, for use on the end. The parameter f controls the parity of the 
function and is zero for even parity, unity for odd parity. Another function, involving the 
polynomial alone, is given by 

^,r,^H ( - ir21+ 'g^^ um-rn, m + V+ y + ?;y + W). (3) 

The multipliers in (2) and (3) are designed to simplify both the orthogonality relation and 
certain integral transforms (see app. A) which play an important part in the analysis. The 
former is simply 

j ^nU 2y - 1 du=8 mn . (4) 

The polynomials themselves are adaptations of Gegenbauer and Jacobi polynomials ; expressions 
for \p m and \f/ m in terms of Jacobi polynomials of argument (2u 2 — 1) are possible. 

3. The Longitudinal Problem 

In the longitudinal problem, a uniform electric field of magnitude E is applied in the 
positive ^-direction. The assumed charge density on the side depends upon the coefficients r b 
and r m , 0<m< N s —1, as follows: 

■V*«- l rb ^ + S r J '"(l h ~if) w**l ( 5) 

€oE I o, |.|>U 

The linear term, 

is called the "basic" term and is included in order that (5) may more easily approximate the 
true distribution especially when a/b is small. 

The potential generated in the interior by the charge density (5) acting alone is easily 
expressed as an integral transform : 



V s (p,z)__a 



cE 



= c f " \l f R ^ sin kz ' dz '~\ KoiMUkp) sin kz dk. (7) 



The Fourier sine transform of the charge density, enclosed in square brackets above, may be 
found with the aid of appendix A. The function I (kp) sin kz, which is regular in a neighbor- 
hood of the origin, may be represented as a spherical harmonic expansion about the origin. 
Thus (7) becomes: 

Vs(.r,e) .^a6/2V" fT r /s#) ^' (- 1)- J„ +>/ , + ,(M) -fc „ . (tcY ^ df . 
cE -£j c W Jo L' W" ,+ ^o" (kb) m+ ' J*ot*on*c; <ut 



(2P+1) 



y P 2p+1 (cos0). (8) 

This expansion converges within that sphere which has its center at the origin and which 
passes through the nearest point occupied by charge, in other words, for all r<a. 

Consider next the ends of the cylinder. The assumed charge density on the upper end 
depends upon the coefficients t b and t m 0<m<N e —l, as follows: 

^T^J^+So^H^'-hO' P - a \ (9) 

I 0, P >a. J 
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The constant term, 

< 6 =^?o(l,0,0,£)> (10) 

is now the "basic" term. Again, the potential generated in the interior of the cylinder by both 
the charge density (9) and an equal and opposite density situated on the lower end, is easily 
expressed as an integral transform : 



V e {p,z)_a 



(" R r T(p') J (kp') P 'dp'~\e- k >>Jo(kp) sinh kz dk. (11) 



This, in turn, may also be given as an expansion in spherical harmonics which converges for 
all r<b: 

cE -p^cjo [_» ka + ^ m (ka) l+ > ] e m dk (2p+iy\c) 'Vh (cosfl). 

(12) 

Evidently the sum of (8) and (12) constitutes an expansion of the interior potential due 
to the total charge distribution. Although this expansion converges only within the largest 
sphere which can be inscribed in the cylinder, it can be continued [3, p. 196] to any interior 
point of the cylinder. This expansion of the interior potential may be expressed in the follow- 
ing form: 

^^P-S nX s J>+ S r m X s p m +t b X p b + S t m XT ( I ) P 2v+l (cos d), (13) 

where the X's are obtained from (8) and (12) by integrating [4, p. 137] over k. These X's 
become the matrix elements in a system of simultaneous linear equations which may be solved 
for the r b , r m , t b , and t m such that the coefficient of (r/c) 2p+1 P 2p +i(cos 0) will be unity for p = 
and zero for as many #>0 as possible. A solution of this type approaches exactness as the 
number of simultaneous equations approaches infinity; however, good accuracy was obtained 
with only eighteen equations. Among these is included an additional relationship known as 
the "edge condition," 



S/i»-(t) S*»=0, (14) 



m=0 \ u / m = 



which insures that the side and end densities match one another as they both tend to infinity 
at the edge. The expressions for the matrix elements are as follows: 



X: 



sm_ 



(_!)*+» 2 2 *-T (m+p+|)r (m+3»+|) ^yjyn+2 

r (|)r (2^+2)r (2m+§+,) ^ c ^ 

•^(m+p +| m-p+l+v; 2m+|+„; ^); (15) 



(-1)«2 2 '-T (m+p+^)T(m+p+l) &N 
T Q) T(2p+2) r (2m+2+ y ) ^ 

■ i F 1 (m+p+^, m-p±l+v; 2m+2+p; J)- (16) 

Separate expressions for Xp and Xf "are not necessary since these elements may be obtained 
merely by setting v=0 in (15) and (16). 
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The dipole moment in the longitudinal problem has only a 2-component and is found by 
integrating za over the entire surface of the cylinder. The orthogonality properties of the tm 
and \p m greatly simplify this integration, which resnlts in the following formula for the longi- 
tudinal polarizability : 



«ee = 26[~l 

» e a \jl 



T(3/2) 



L 3^+2 2/:i r(13/6) ro J 



'°J +2 U** + 2 2 ' 3 r(5/3)* J' 



(17) 



where v is the geometrical volume of the cylinder. 

4. The Transverse Problem 

In the transverse problem, a uniform electric field of magnitude E is applied in the positive 
^-direction. The assumed charge density on the side depends upon the coefficients s b and s m , 
0<m<N s —l, as follows: 



e E 



= [cos </>]£(£) = [cos 0]-< 



(18) 



form 



M-l - /I 1 S\ 

J), \*\>b. m 

The potential of tins density in the interior of the cylinder is given by the integral trans- 

^^=| r R fsiz') cos kz' dz'~\ K^kailxikp) cos kz dk cos*. (19) 

Its spherical harmonic representation, convergent for r<^a, becomes 
V.irM) ^ab/2\» f T Jmikb) N ^\ (-l)»J im+m+ ,W) l 

'Pic W Jo L* (kby /i+ i=i> m m m+ " J 



cE 



cE 



■KMaKkcY'+idk^j^^* 1 P\ p +1 (cos 6) cos*. (20) 

The assumed cliarge density on the upper end involves the coefficients w b and w m , 0<m 

<N e —l, and is given by: 



€0 



E 



= [cos 4>W(p) = I cos 0] * 



p Ne-l _ / I p \ 

Wr+SW 1, 1, —7V-)>P<<L. 

a m =o \ 6 a/ 

, p>a. 



(21) 



The potential generated in the interior of the cylinder by both this charge density and an equal 
densitv on the lower end becomes 



V e (p,<t>,z) 
cE 



--- PR p W( P ') JiCipOp'rfp'l e-^Mkp) cosh kzdk cos <£. (22) 



Again, the spherical harmonic expansion, convergent for r<6, is given by 
V e (r,d,<p) _^\a 2 f«T J2W , V ,„ (-l) w c/2m + 2 + >(to) "l 

c# ~^ocJ L 6 to ^^i m ~ (*a) 1+F J 



1 / r \2p+l 

e -m (ke y P+ i dk _^_^ pi p+1 (co3 9)cos*. (23) 
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Once more, the sum of (20) and (23) constitutes an expansion of the interior potential 
due to the total charge distribution, which expansion converges within the largest sphere that 
can be inscribed in the cylinder. The expansion may be expressed as follows: 



V(r 9 e,4>)_ 



cE 



oo r Ns-l Ne-l H / r \2p + l 

=1] s b Y s p b + S 8 n Y' p m +w b Y?+ 2 w m Y™ ( - ) iVi(cos 9) cos 0, (24) 

p=0 L ra=0 ra=0 J \C/ 



where the F's are again obtained by integrating over k. It is desired that s b , s m , w b , and w m 
be such that the coefficient of (r/c) 2p+1 P 2p +i(cos 6) cos # be unity for p = and zero for as many 
p^>0 as possible, and a system of simultaneous linear equations must be solved. The edge 
condition now becomes 



N '~ 1 //7\ 1/3 Ne-l 

w=0 \0/ w=0 



and the expressions for the matrix elements are: 



w„ 



=0, 



(25) 



ysm_ 

J- p - 



(-D"-a*-T (m+i»+|) r ( m+i ,+i) 

WW 



J\2w+2 



r(i)r(2p+3)r(2m+|+ P ) 



/ 3 3 5 2 \ 



(26) 






rQ)r(2^+3)r(2m+3+^) 



(27) 



Note that Y v b and Yp b are obtained by setting *>=0 in the above. 

The dipole moment now has only an x component and is found by integrating xa over 
the entire surface of the cylinder. The expression for the transverse polarizability is: 






r(i/2) 



;« 



1 



4 W;&n "2 2/3 r(8/3) 



W ' 



(28) 



2 2/3 r(7/6) 

5. Results 
The two distinct tensor components of the electric polarizability are given in table 1. 

Table 1 



a/6 


an/floeo 


att/voeo 





00 


2. 0000 


H 


15. 071 


2. 3151 


H 


7.0966 


2. 6115 


l 


3. 8614 


3. 1707 


2 


2. 4325 


4. 2173 


4 


1. 7507 


6. 1814 


CO 


1. 0000 


CO 



VQ=2Tra*b 



These results may be compared with the polarizability of a conducting sphere which is equal 
to 3v e in rationalized M.K.S. units. Actual values of the charge coefficients are in- 
cluded in appendix B. An estimate of the accurac}^ of the results was obtained by solving 
the problem repeatedly with larger and larger numbers of equations and observing the limit 
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toward which the calculated value of the polarizability tended. It was found that, for a maxi- 
mum number of equations equal to 18, the polarizability could he found to 5 significant figures 
with some doubt about the least significant figure. As might be expected, the accuracy is best 
for a/b equal to unity, in which case the least significant figure is in error by at most one unit. 
As an additional check upon the calculations, the potential, V, at the pole and the electric 
field, E z , at the equator due to the charge alone in the longitudinal problem were calculated 
and compared with the corresponding values associated with the applied field. The small 
errors involved can, by appealing to the appropriate gradients, be translated into local dis- 
crepancies, Ap at the equator and Az at the pole, between the actual surface and an imaginary 
surface which is an equipotential under the joint influence of charge and applied field. A 
corresponding check, in which the potential, V, at the equator and the electric field, E p , at the 
pole due to the charge alone were compared with the corresponding values associated with the 
applied field, was carried out for the transverse problem. The results of both checks are given 
in appendix B and are seen to be very satisfactory. 

6. Appendix A: Hankel Transforms of the \p m Functions 

It will be shown that all of the integral transforms of the \(/ m functions used in this article 
are special cases of the general transform: 



J. 



V(l-^)^>(2u^ 



Re o-> — 1; Re *>>-!; t real and >0, (Al) 



which will now be proved. P^' a) (2u 2 —l) is a Jacoby polynomial [5, p. 168]. Denote the in- 
tegral in eq (Al) hyf(t). Then 

Letting 2u 2 — l=x in [4, p. 284], 



x 



\>XX-uyP^{2u*-l)udu= r i^^ 7 r(m+1 + ^T" +1) , , ;Re*>-l;Re,»-l. 



2r(m+l)r(J— <r— m+l)T(v+l+m+2) 
Let l=r-\-(r; Re o-> — 1: 



(A3) 



r(m+l+„) ^ (-l) f (f/2)»+> „ 

J(t)= , rrrr- > , — , iNt ./ — : ; ; r- ^r, Ke ff> — 1 , K(5 V>— 1. (A4) 

2r(m+l) ^cr(r-m-(-l)r(i'+r-fff+m+2)' 
Make a change in the index of summation; let r=m-\-k: 

f m = (- 1 ) mT ( m + 1 + v ) ^ (-i)*(f/2)*+»»+* 

]W 2r(m+i) £ir(i+ic)r(2m+2+a+ v +ky K o; 



This is: 



f{t) 2r(m+l) W «W»+'+'W> ( A6 ) 



and the transform of eq (Al) is true. 
From [5, p. 170], one finds that: 



P<-> (2^-l) = ( T ^^^ 2 F 1 (-m,m+v+,+ l;<r+V,u^ (A7) 
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If o is set equal to 7 + f— 1, the transform becomes 



£ Kh,S,v,v) Jy+t-Mu) {tu)Hu= ( 1)m f ^-y+J + 



ay 



(A8) 



Now if the parameters 7 and f are allowed to take all values relevant to the present problem, 
one obtains: 



Side, even: 



Side, odd: 



End, even: 



End, odd: 



o f m \2'°> V > M ) cos tudu =\^) 
j tm \^\ } v,u\ sin tu du=(^\ - — 



r\ 1,a (-l) W e7 2w+ i/2 + > (Q . 

^1/2+y 



( l) m J 2m±l + v\t) . 



£ MhO,»,u)J (tu)udu= ( ir ^Z 
£jMA,r',n)J 1 (tu)udu^ ( — ir fz +2+y{t) > ^ v>-\. 



(A9) 
(A10) 
(All) 
(A12) 



Because of the special values taken by 7 and f , the singularity at the origin which might have 
occurred in (tu) 7 J y+ ^i(tu) does not appear. Both sides of eqs (A9) through (A12) become 
entire functions of t and no restrictions need be placed upon t. 

7. Appendix B: Values of the Charge Coefficients 

The values of the charge coefficients and the results of the checking procedure at pole and 
equator are given in tables 2 through 6. The notation "A(p)" means "A times lO 2 * ,; 

Table 2 Table 3 



a/b = U 





Longitudinal 


Transverse 


m 


r m 


Sm 


Basic 



1 
2 
3 
4 

5 
6 

7 
8 
9 

10 
11 
12 
13 


+0.12905187 (+0) 

+0.23603290 (+1) 
-0.11393541 (+0) 
-0.19185823 (-1) 
-0.52174681 (-2) 
-0.16484165 (-2) 

-0.53665621 (-3) 
-0.16813125 (-3) 
-0.48283874 (-4) 
-0.12195426 (-4) 
-0.25981907 (-5) 

-0.44409406 (-6) 
-0.56727317 (-7) 
-0.47918712 (-8) 
-0.20011235 (-9) 


+0.40387637 (+0) 

+0.14690787 (+1) 
-0.43239914 (+0) 
-0.65706908 (-1) 
-0.17221441 (-1) 
-0.54219231 (-2) 

-0.17807221 (-2) 
-0.56407302 (-3) 
-0.16353609 (-3) 
-0.41591873 (-4) 
-0.88990999 (-5) 

-0.15241129 (-5) 
-0.19470442 (-6) 
-0.16423381 (-7) 
-0.68402739 (-9) 


Equator 
check 


EzlE= -1.00000002 


V/aE cos <f>= +1.00000002 


Ap/a 


Negligible 


Negligible 


m 


tm 


w m 


Basic 


1 


+0.15200248 (-1) 

+0.35624820 (+1) 
-0.39128623 (-1) 


-0.17172000 (+0) 

+0.14693027 (+1) 
+0.32011553 (-1) 


Pole check 


^lbE= +0.99999875 


Ep/E cos <f>= -0. 99902328 


Az/b 


-0.00000037 


-0. 00022077 



a/b=M 





Longitudinal 


Transverse 


m 


r m 


Sm 


Basic 


1 
2 
3 
4 

5 
6 

7 
8 
9 
10 


-0.19225150 (-4) 

+0.18367917 (+1) 
-0.48838691 (-1) 
-0.54930270 (-2) 
-0.10582670 (-2) 
-0.23834773 (-3) 

-0.54147255 (-4) 
-0.11308514 (-4) 
-0.20056279 (-5) 
-0.27611488 (-6) 
-0.25757920 (-7) 
-0.12042175 (-8) 


-4-0.38965487 (+0) 

+0.16055212 (+1) 
-0.35509127 (+0) 
-0.43967119 (-1) 
-0.10119125 (-1) 
-0.28016839 (-2) 

-0.78539136 (-3) 
-0.20187101 (-3) 
-0.43956298 (-4) 
-0.74296836 (-5) 
-0.85391338 (-6) 
-0.49561337 (-7) 


Equator 
check 


#,/#=- 1.00000003 


T '/aE cos </>=+l. 00000009 


A P /a 


Negligible 


Negligible 


m 


tm 


Wm 


Basic 


1 
2 
3 
4 


+0.46938430 (-1) 

+0.23357609 (+1) 
-0.82663119 (-1) 
-0.76227915 (-2) 
-0.12692924 (-2) 
-0.16584913 (-3) 


+0.50835705 (-2) 

+0.13797193 (+1) 
+0.88406931 (-1) 
+0.25064651 (-1) 
+0.77264312 (-2) 
+0.15415744 (-2) 


Pole check 


^JbE= +0.99999997 


Ep/E cos 0=-l. 00002004 


A zjb 


Negligible 




+0. 00000890 
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T \BLE 4 



a/b = l 





Longitudinal 


Transverse 


m 


Tm 




Sm 


Basic 


1 
2 
3 

4 
5 
6 

7 


-0.37616817 (-1) 

+0.14921036 (+1) 
-0.83747159 (-2) 
+0.23944987 (-2) 
+0.12471808 (-2) 

+0.48160777 (-3) 
+0.14698364 (-3) 
+0.31457070 (-4) 
+0.34949152 (-5) 


+0.25340314 (+0) 

+0.19371617 (+1) 
-0.32219422 (+0) 
-0.38977600 (-1) 
-0.93990090 (-2) 

-0.25945051 (-2) 
-0.64974220 (-3) 
-0.12162489 (-3) 
-0.12244613 (-4) 


Equator 
check 


E z /E= -1.00000001 


* r /aE cos = +0. 99999999 


Apia 


Negligible 


Negligible 


m 


tm 


w m 


Basic 


1 
2 
3 

4 

5 
6 

7 


+0.87113326 (-1) 

+0.16193448 (+1) 
-0.11224086 (+0) 
-0.14340444 (-1) 
-0.34812949 (-2) 

-0.95870761 (-3) 
-0.23986408 (-3) 
-0.44962281 (-4) 
-0.45407822 (-5) 


-0.13425511 (+0) 

+0.14792105 (+1) 
+0.58149961 (-1) 
+0.17095607 (-1) 
+0.59774893 (-2) 

+0.20469643 (-2) 
+0.59421158 (-3) 
+0.12429137 (-3) 
+0.13669068 (-4) 


Pole check 


V/bE= +1.00000001 


Ep/E cos </>= -0. 99999996 


A 2/6 


Negligible 




Negligible 





Table 5 



ajb =2 





Longitudinal 


Transverse 


m 


r m 


Sm 


Basic 


1 
2 
3 
4 


-0.65161515 (-2) 

+0.12482169 (+1) 
+0.24818534 (-1) 
+0.91060187 (-2) 
+0.28707901 (-2) 
+0.56748900 (-3) 


+0.16191166 (+0) 

+0.23744087 (+1) 
-0.28135556 (+0) 
-0.33694020 (-1) 
-0.72222609 (-2) 
-0.11581174 (-2) 


Equator 
check 


E t jE= -1.00000271 


V/aE cos «/» = +0. 99999988 


Apia 


+0. 00000121 


-0. 00000005 


m 


U 


w m 


Basic 


1 
2 
3 
4 

5 
6 
7 
8 
9 
10 


+0.16025169 (+0) 

+0.11683321 (+1) 
-0.12492981 (+0) 
-0.17150631 (-1) 
-0.42141901 (-2) 
-0.12098644 (-2) 

-0.34737952 (-3) 
-0.90947390 (-4) 
-0.20111927 (-4) 
-0.34459642 (-5) 
-0.40094338 (-6) 
-0.23533528 (-7) 


-0.22721442 (+0) 

+0.15658591 (+1) 
+0.42693596 (-1) 
+0.12814464 (-1) 
+0.43243348 (-2) 
+0.14993639 (-2) 

+0.49000507 (-3) 
+0.14126899 (-3) 
+C. 33697692 (-4) 
+0.61418057 (-5) 
+0.75263582 (-6) 
+0.46182535 (-7) 


Pole check 


Vfb E= +1.00000002 


Ep/E cos <t> = -1 . 00000003 


Azjb 


Negligible 


Negligible 



Table 6 



a/b=i 





Longitudinal 


Transverse 


III 


Ttn 


Sm 


Basic 



1 


-0.91407486 (-1) 

+.11655470 (+1) 
+.92149445 (-2) 


+0.36981774 (+0) 

+.27748758 (+1) 
-.13470990 (+0) 


Equator 
check 


Ez/E= -0.99979406 


VlaE cos 0=+l. 00001613 


Apia 


-0.00005287 


+0.00000530 


m 


U 


w m 


Basic 


1 
2 
3 
4 

5 



7 
8 
9 

10 
11 
12 
13 


+0.16831628 (+0) 

+0.92928755 (+0) 
-0.15311162 (+0) 
-0.25352751 (-1) 
-0.71931540 (-2) 
-0.23904452 (-2) 

-0.81543454 (-3) 
-0.26584520 (-3) 
-0.78893956 (-4) 
-0.20466149 (-4) 
-0.44553523 (-5) 

-0.77488570 (-6) 
-0.10037598 (-6) 
-0.85747540 (-8) 
-0.36132625 (-9) 


-0.10545599 (+0) 

+0.15388757 (+1) 
+0.82957625 (-1) 
+0.25448502 (-1) 
+0.96857682 (-2) 
+0.39007948 (-2) 

+0.15309282 (-2) 
+0.55035828 (-3) 
+0.18020312 (-3) 
+0.50262768 (-4) 
+0.11635018 (-4) 

+0.21335661 (-5) 
+0.28945170 (-6) 
+0.25757752 (-7) 
+0.11250670 (-8) 


Pole check 


176E=+1.00000C01 


Ep/E cos <t> = -1.00000001 


Azjb 


Negligible 


Negligible 
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